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Abstract. We have adopted the multistep shell model in the complex energy plane to study 
nuclear excitations occurring in the continuum part of the spectrum. In this method one 
proceeds by solving the shell model equations in a successive manner. That is, in each step 
one constructs the building blocks to be used in future steps. We applied this formalism to 
analyze the unbound nuclei 12,13 Li starting from the one-particle states in 10 Li and two-particle 
states in 11 Li. In the former case the excitations correspond to the motion of three particles 
partitioned as the product of a one-particle and two-particle systems. The ground state of 12 Li 
is thus calculated to be an antibound (virtual) state. In the four-particle system 13 Li the states 
can be constructed as the coupling of two correlated pairs. We found that there is no bound or 
antibound state in 13 Li. 



1. Introduction 

Experimental facilities allow one nowadays to measure systems living a very short time. To 
describe these processes one has to consider the decaying character of the system. On the other 
hand, the theoretical study of unstable nuclei is a difficult undertaking. However, the system 
may be considered stationary if it lives a long time. In this case the time dependence can be 
circumvented. Of the various theories that have been conceived to analyze unbound systems, we 
mention an extension of the shell model to the complex energy plane pQ . The basic assumption 
of this theory is that resonances can be described in terms of states lying in the complex energy 
plane. The real parts of the corresponding energies are the positions of the resonances while 
the imaginary parts are minus twice the corresponding widths, as it was proposed by Gamow 
at the beginning of quantum mechanics [2]. These complex states correspond to solutions of 
the Schrodinger equation with outgoing boundary conditions. Details of the formalism and its 
application to two-nucleon systems can be found, e.g., in Refs. 0H1E]. 

Our aim is to develop a suitable formalism to treat unstable nuclei involving many valence 
nucleons in the continuum [6 J . This formalism is an extension of the shell model in the complex 
energy plane [T]. The correlations induced by the pairing force acting upon particles moving in 
decaying single-particle states is taken into account by using the multistep shell model (MSM) [7j. 
The formalism is presented in Section [2J Applications are in Section and a summary and 
conclusions are in Section 01 
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Figure 1. Integration contour L + in the complex energy plane. The open circles denote the 
resonances included in the sum of Eq. 0, while the solid circles are those excluded. The vertex 
(c, 0) corresponds to the energy cutoff point c. 



2. The formalism 

The eigenstates of a central potential obtained as outgoing solutions of the Schrodinger equation 
can be used to express the Dirac <5-function as [8] 

*(r-0 = 5>M*.M+/ (i) 

n J L+ 

where the sum runs over all the bound and antibound states plus the complex states (resonances) 
which lie between the real energy axis and the integration contour L + (see, Fig. [1]). The 
antibound states are virtual states with negative scattering length. The wave function of a state 
n in these discrete set is w n {r) and u(r,E) is the scattering function at energy E. 

Discretizing the integral of Eq. ([1]) one obtains the set of orthonormal vectors \ipj) forming 
the Berggren representation [9]. Since this discretization provides an approximate value of the 
integral, the Berggren vectors fulfill the relation I ~ • \(fj){<fj\- One has, 

/ dEu(r, E)u(r', E) = V h p u(r, E p )u(r' , E p ), (2) 
Jl+ p 

where the quantities E p and h p are defined by the procedure that uses to perform the 
integration. Therefore the orthonormal vectors \(fj) are given by the set of bound, antibound 
and Gamow states, i.e., (r\(p n ) = {w n (r, E n )}, and the discretized scattering states, i.e., 



(r\<p p ) = {y/h^,Up(r, E p )}. In Ref. [9] it was also found that few discretized scattering states in 
the basis are enough to obtain convergence. Using the Berggren representation one readily gets 
the two-particle shell-model equations in the complex energy plane (CXSM) [3], i.e., 

(W(a 2 ) -€i- ej)X(ij; a 2 ) = J^kl; a 2 \V\ij; a 2 )X(kl; a 2 ), (3) 

k<i 

where V is the residual interaction and i,j,k,l label single-particle states. The two-particle 
states are labeled by a 2 . e and W denote the energy of the single-particle and two-particle 
states, respectively. The tilde in the interaction matrix element denotes mirror states so that in 
the corresponding radial integral there is not any complex conjugate, as required by the Berggren 
metric. 

The two-particle wave function is given by 

\a 2 )=P+(a 2 )\0), (4) 
where the two-particle creation operator is 

P+(a 2 ) = J2x(ir,a 2 ) iC )f^ , (5) 

and X a2 is the angular momentum of the two-particle state. 
2.1. Separable interaction 

For a separable interaction the matrix can be written as 

(kl;a\V\ij;a) = -G a f a (kl)f a (ij), (6) 

where G a is the interaction strength and f a (ij) is the matrix element of the field defining the 
interaction. The shell model equation becomes 

_ X(ij;a) = f a (kl)f a (ij)X(kl;a) 

G a ~ {u a -a- ej) 

Multiplying Yli<j fa(ij) on both sides of Eq. ([7]) one gets the so-called dispersion relation 



<Jq . ~ . U a — 6j — €j 



The two-particle wave function amplitudes are given by 



X(ij;a)=N a Uij) , (9) 

(jJ n — €i — €-i 



where N a is the normalization constant determined by requiring 



^2x\ij;a) = l. (10) 



2.2. The Multistep Shell Model Method 

The Multistep Shell Model Method (MSM) solves the shell model equations in several steps. In 
the first step the single-particle representation is chosen. In the second step the energies and 
wave functions of the two-particle system are evaluated with a given two-particle interaction. 
The three-particle states are evaluated in terms of a basis consisting of the tensorial product 
of the one- and two-particle states previously obtained. The MSM basis is overcomplete and 
non-orthogonal. To correct this one needs to evaluate the overlap matrix among the basis states 
also. A general description of the formalism is in Ref. [TU]. The particular system that is of 
our interest here, i.e., the three-particle case, can be found in Ref. [7]. Below we refer to this 
formalism as CXMSM [BJ. 

The three-particle energies VF(a3) are given by [7] 

(W(a 3 ) -Si- W{a 2 )){a 3 \(cfP+{a 2 )) a3 \0) 
= E I E(^) -ei- e k )A(ia 2 ,jp 2 ; k) \ (a 3 \(c+P + (p 2 )) a3 10), (11) 

where 

A(ia 2 , j p 2 ;k) = a 2 p 2 {] ^ M Y(kj;a 2 )Y(ki;j3 2 ), (12) 

and 

Y(ij;a 2 ) = {l + 5{i,j)) l l 2 X{ij;a 2 ). (13) 

The rest of the notation is standard. 

The matrix defined in Eq. (jlip is not hermitian and the dimension may be larger than the 
corresponding shell-model dimension. This is due to the violations of the Pauli principle as well 
as over-counting of states in the CXMSM basis. Therefore the direct diagonalization of Eq. (jlip 
is not convenient. One needs to calculate the overlap matrix in order to transform the CXMSM 
basis into an orthonormal set. In this three-particle case the overlap matrix is 

(0\(ctP + (a 2 ))l i (cjP + (f3 2 )) as \0) =8 lJ 8 a2 p 2 + y £A(ia 2 ,j(3 2 ;k). (14) 

k 

Using this matrix (|14p one can transform the matrix determined by Eq. (jlip into a hermitian 
matrix T which has the right dimension. The diagonalization of T provides the three-particle 
energies. The corresponding wave function amplitudes can be readily evaluated to obtain 

\a 3 ) = P + (a 3 )\0), (15) 
P + (a 3 ) = ^X(fa 2 ;a 3 )(c+P + (a 2 )) Q 3, (16) 

ia.2 

where P + {a 3 ) is the three-particle creation operator. It has to be pointed out that in cases 
where the basis is overcomplete the amplitudes X are not well defined. But this is no hinder to 
evaluate the physical quantities. For details see Ref. [7]. 

The CXMSM allows one to choose in the basis states a limited number of excitations. This 
is because in the continuum the vast majority of basis states consists of scattering functions. 
That is, the majority of the two-particle states provided by the CXSM are complex states which 
form a part of the continuum background. Only a few of those calculated states correspond 
to physically meaningful resonances, i.e., resonances which can be observed. Below we call a 
"resonance" only to a complex state which is meaningful. These resonances are mainly built 
upon single-particle states which are either bound or narrow resonances. Yet, one cannot ignore 



the continuum when evaluating the resonances. The continuum configurations in the resonance 
wave function are small but many, and they affect the two-particle resonance significantly [1]. 
Therefore, the great advantage of the CXMSM is that one can include in the basis only two- 
particle resonances, while neglecting the background continuum states, which form the vast 
majority of complex two-particle states. 

2.3. The neutron-proton correlation 

The influence of the neutron-proton correlation can be considered in a straightforward way in 
our CXMSM. For example, we consider a system with three neutrons and one proton , say 12 Li 
by assuming 8 He as the inert core [6]. With the basis denoted as 

\pa 3 ;a 4 ) = (c+P+(a 3 )) a4 10), (17) 

where p labels the proton state and are the three-neutron one-proton state, the four-particle 
energies VK(«4) in 12 Li are given by 

(W(a 4 ) -e p - W(a 3 )) (a 4 |(c+ P+(a 3 )) Q4 |0) 
= E| E^ A I^ A > Bi+ E ^';A|y|^;A)B 2 l(«4|(c+P + (^))a4|0), (18) 

<j/?3 yklXa-2 ijkl\a 2 /32 ) 

where, 

a = |-ir + '^(»« S «,W«.;AfeW ! { £ * Q A s }{ I J 4 I }, (19) 

and 

B 2 = (-l)P + ^ i+l Y(ji;a 2 )Y(jk;P 2 )X(ka 2 ;a 3 )F(W^ 3 ) 

xa 2 a 3 /y 3 A 2 { I 1 A U * 1 X R ){ I j " 2 ). (20) 
{ p 2 a 4 a 3 } { a 4 p 3 ) \ k a 3 j3 2 J 

Here p and q label proton states, while i,j, k, I label neutron states, (pk; \\V\ql; A) denotes the 
corresponding proton-neutron interaction matrix element. The wave function amplitudes X and 
the projected quantities F have been evaluated in previous steps of the CXMSM. Notice that 
in this case the overlap matrix is the unit matrix, i.e., 

(0\(c+P+(a' 3 ))i A (c+P + (a 3 )) ai \0) = V<W (21) 

2. 4- Four-particle states 

Systems with four like particles can be described in a correlated two-particle basis as 

|« 2 /3 2 ;a 4 ) = (P t (a 2 )P+(/3 2 )) Q4 |0), (22) 
where a 2 ^ f3 2 . One can obtain the four-particle MSM equation as 

(W(a 4 )-W(a 2 )-W(P 2 ))(a 4 \(PHa 2 )PH^))JO) = - £ £(1 + A^)" 1 

72<<52 ijkl 

x(W( 72 ) + W{5 2 ) -a- €j -e k - ei)B(ijkl, a 2 p 2l2 5 2 ; a A ){a A \(P\ l2 )P\5 2 )) aA \U) , (23) 



where 



i j a 2 

B(ijkl,a 2 p 2 j2S2]u 4 ) = a 2 M2S2Y(ij;a 2 )Y(kl-,p 2 )Y(ih, l2 )Y(jl-,5 2 ) { k I fa } . (24) 

72 5 2 04 

The four-particle overlap matrix can be obtained in a similar way as 

(0|(pt( a2 )pt (/32)) t 4( pt (72)P t ((52))a j 0) = Sa2l2 6p 252 

+ (-lT 2+ ^- a4 Sa 2 S 2 Sp 2J2 - Yl B ^ kl ' ^2l2h). (25) 

ijkl 

To restore the Pauli principle, one can use the Schmidt procedure to generate a set of 
orthonormal basis \ui) (which has the right dimensions) from the overlap matrix. The overlap 
matrix is thus transformed into an identity matrix 

(Um\u n ) = Iran = ^ €m(i)0(i, j)£n(j) > ( 26 ) 
ij 

where 0(i,j) is the overlap matrix, and the wave function in the Schmidt basis is 

\ u m) = /^m(i)\i}- (27) 

i 

One can transform the dynamical matrix into a hermitian matrix T in the right dimensions if 
all possible basis states are included. The dynamical matrix M(ij) in the MSM basis is given 
by 

W(a)(i\a) =Y f M(i,j){j\a). (28) 
j 



To transform it into the orthonormal basis, one can use Eq. (|27p and the identity equation. One 
thus obtains 

W(a)Y,Hm(i)(i\a) = ^Cm« £ M(i, j) ^(j|n n )(u„|a), (29) 

i i j n 

W(a)(u m \a) 

n ^ ij k * 

The expression in the bracket is therefore the hermitian matrix T, which can be written as 

T(m,n) = Y,U(i)M(i,j)0(j,k)Uk). (31) 

ijk 



3. Applications 

Using the Berggren single-particle representation described above, we will evaluate the complex 
energies and wave functions of the unbound nuclei 12,13 Li using the MSM basis states consisting 
of the Berggren one-particular states, which are states in 10 Li, times the two-particle states 
corresponding to 11 Li. The spectrum of n Li was already evaluated within the CXSM including 
antibound states |13j . Here we will repeat that calculation in order to determine the two-particle 
states to be used in the calculation of the three- and four-particle systems. 

To define the Berggren single-particle representation we still have to choose the integration 
contour L + (see Eq. {!])). The valence shells are the low lying resonances Opi/2 at (0.195,-0.047) 



MeV and 0d 5 / 2 at (2.731, -0.545) MeV. Besides, the state lsi/ 2 appears as an antibound state. 
To include in the representation the antibound lsi/2 state as well as the Gamow resonances Opi/2 
and Ods/2 we will use two different contours. The number of points on each contour defines the 
energies of the scattering functions in the Berggren representation, i.e., the number of basis 
states corresponding to the continuum background. This number is not uniformity distributed, 
since in segments of the contour which are close to the antibound state or to a resonance the 
scattering functions increase strongly. 
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Figure 2. Left: Contour used to include the antibound state (see, also, Ref. [T2] ). The points 
Bi denote bound state energies while A denotes the antiboud state. Right: Contour used to 
include the Gamow resonances represented by the points Gj. 



The dynamics of Li is determined by the pairing force acting upon the two neutrons coupled 
to a state + , which behaves as a normal even-even ground state [IT] [15] . Besides the energy, 
this state has been measured to have an angular momentum contain of about 60 % of s-waves 
and 40 % of p-waves, although small components of other angular momenta are not excluded. 

We will use a separable interaction as in Ref. [13]. The strength G\ 2 , corresponding to 
the states with angular momentum A2 and parity (— 1) A2 , will be determined by fitting the 
experimental energy of the lowest of these states, as usual. It is worthwhile to point out that 
G\ 2 defines the Hamiltonian and, therefore, is a real quantity. The energies are thus obtained 
by solving the corresponding dispersion relation. The two-particle wave function amplitudes are 
given by 

Yf ^ AT /fa'. «2) / Q0 x 

X(ij;a 2 ) = N a2 — ■ -, (32) 

where f(ij, a 2 ) is the single particle matrix element of the field defining the separable interaction 
and N a2 is the normalization constant determined by the condition Y2i<j a 2 ) 2 = 1. 

Due to the large number of scattering states included in the single-particle representation 
the dimension of three-particle basis is also large. The scattering states are needed in order 
to describe these unstable states. In the calculations we took into account all the possibilities 
described above regarding the energies of the single-particle state 0pi/ 2 as wen as the binding 
energy of the state 11 Li(gs). 

With the single-particle states and the two-particle states 11 Li(gs) and n Li(2f) discussed 
above, we formed all the possible three-particle basis states. We found that the only physically 
relevant states are those which are mainly determined by the bound state n Li(0^). The 



corresponding spins and parities are l/2 + , 1/2 and 5/2 + . States like 3/2 + , which arises from 
the CXMSM configuration |ls 1/2 ® 2^;3/2 + ), is not a meaningful state. The corresponding 
calculated energies are shown in Fig. [3j 




Figure 3. Experimental level scheme in 12 Li. The three lowest levels are from [16], while the 
one at 1.5 MeV is from [15]. In the second column are the three-neutron CXMSM results. In 
the columns A-C are the shell-model calculations corresponding to different truncation schemes: 
A) Maximum of lp excitation from p to sd shell, B) maximum of 3p excitations and C) full psd 
space. Dashed lines indicate the widths of the resonances. 

In Ref. [16] it was also found that 12 Li(gs) is an antibound state but, in addition, two 
other low-lying states were observed at 0.250 MeV and 0.555 MeV by using two-proton removal 
reactions. In this case the 0^3/2 proton in the core may interfere with the neutron excitations 
evaluated above. In particular, the antibound 1/2+ ground state would provide, through the 
proton excitation, a state 1~ and a 2~. This is the situation encountered in the shell model 
calculation [17]. One sees in this Figure that the full calculation predicts all excited states to 
lie well above the corresponding experimental values. It is worthwhile to point out that the 
calculated states exhibit rather pure shell model configurations. For instance the states 2^~ 
(ground state) and l± are mainly composed of the configuration \n [0p 3 / 2 ] v [(0p 1 / 2 ) 2 l'Si/2] )• 
This does not fully agree with our CXMSM calculation, since in our case this wave function is 
mainly of the form |lsi/2 <S> n Li(gs)). This differs from the shell model case in two ways. First, 
the state 11 Li(gs) contains nearly as much of lsi/2 as of 0p\/2- Second the continuum states 
contribute much in the building up of the antibound 12 Li(gs) wave function, as discussed above. 
In our representation it is straightforward to discern the antibound character of this state, which 
is not the case when using harmonic oscillator bases. 

Using the two-particle states of n Li within the Berggren single-particle representation, we 
can calculate the four-particle system 13 Li. The four-particle CXMSM basis is constructed by 
two correlated two-particle states. The two-particle states we included here are the 11 Li(0^ 1 "), 
11 Li(l^~) and 11 Li(2^), which are the same as the three-particle case. In this basis set we formed 



Table 1. Calculated energies (in MeV) of the four-particle states in 13 Li corresponding 
to different single-particle e Pl/2 energies. The two-particle states W(0f) = — 0.295MeV, 
W(ir/) = (0.084, -0.002) and W{2f ) = (2.300, -0.372)MeV are included in the calculation. 





e Pl/2 = (0.195,-0.047) 


e Pl/2 = (0.563, -0.252) 


0+ 


(0.868,-0.059) 


(1.505,-0.041) 




(5.127,-0.964) 


(5.274,-1.009) 


1- 


(0.836,-0.117) 


(1.722,-0.166) 




(3.244,-0.593) 


(3.921,-0.721) 


2+ 


(0.715,-0.114) 


(1.802,-0.257) 




(2.907,-0.445) 


(3.373,-0.478) 




(5.131,-0.910) 


(5.205,-0.966) 


3~ 


(2.541,-0.391) 


(2.674,-0.413) 


4+ 


(5.715,-1.119) 


(5.715,-1.119) 



Table 2. Same as Table Q] but for W(0 + ) = -0.369MeV. 





e Pl/2 = (0.195,-0.047) 


e Pl/2 = (0.563, -0.252) 


0+ 


(0.855,-0.057) 


(1.527,-0.033) 




(5.142,-0.970) 


(5.303,-1.020) 


1- 


(0.828,-0.122) 


(1.744,-0.169) 




(3.244,-0.593) 


(3.921,-0.721) 


2+ 


(0.715,-0.114) 


(1.802,-0.257) 




(2.903,-0.450) 


(3.367,-0.477) 




(5.137,-0.916) 


(5.213,-0.973) 


3~ 


(2.541,-0.391) 


(2.674,-0.413) 


4+ 


(5.715,-1.119) 


(5.715,-1.119) 



the symmetric Hamiltonian matrix by evaluating the dynamical matrix Eq. (|23p and the overlap 
Eq. ®. 

The only four-particle states we can get are + , 1~, 2 + , 3~ and 4 + states. The calculated 
results refer to different single-particle e Pl/2 and two-particle 11 Li(gs) energies are all listed in 
Tables [T] and [2j We found that there is no bound or antibound state in 13 Li, which agrees with 
the experiment that the 13 Li is unbound [14]. The ground state is a resonance at around 1 
MeV, however the configuration and the energy strongly depend on which single-particle energy 
e Pi/2 we choose. For e Pl/2 = 0.195MeV case, the ground state is | Li(17/)(g> Li(17/);2 + ) at 
the energy of (0.715, — 0.114)MeV, while for e Pl/2 = 0.563MeV case, the ground state is + 
which at the energy of (1.505, -0.041)MeV and with the components | 11 Li(17;)0 11 Li(l ] ;); 0+) 
and | 11 Li(0^)(8' 11 Li(0^); + ) about half and half. Comparing between tables one can find that 
the four-particle states are not much affected by the energy of 11 Li(gs). 

In our calculation we only included the resonances of 11 Li which is considered physically 
meaningful for building the whole four-particle basis. However as we mentioned before, there 
are lots of other continuum states which belong to the background and have a large component 
of the poles. Although these states are not in themselves physically meaningful resonances, they 
might have some influence on the spectrum of 13 Li. Therefore for a more precise calculation, 
one should include some of these two-particle states. 



4. Summary and conclusions 

We have studied excitations occurring in the continuum part of the nuclear spectrum which are 
at the limit of what can be observed within present experimental facilities. These states are 
very unstable but yet live a time long enough to be amenable to be treated within stationary 
formalisms. We have thus adopted the shell model in the complex energy plane for this purpose. 
In addition we performed the shell model calculation by using the multistep shell model. In 
this method of solving the shell model equations one proceeds in several steps. In each step one 
constructs building blocks to be used in future steps [TO]- We applied this formalism to analyze 
12 ' 13 Li as determined by the neutron degrees of freedom. 

By using single-particle energies (i.e., states in 10 Li) as provided by experimental data when 
available or as provided by our calculation, we found that the only physically meaningful two- 
particle states are 11 Li(gs), which is a bound state, and 11 Li(2^ t ~), which is a resonance. As a 
result there are only three physically meaningful states in 12 Li which, besides the antibound 
ground state, it is predicted that there is a resonance l/2~ lying at about 1 MeV and about 800 
keV wide and another resonance which is 5/2 + lying at about 1.1 MeV and 500 keV wide. That 
the ground state is an antibound (or virtual) state was confirmed by a number of experiments 
|14 |. I15 ^ [T6] and the state 5/2 + has probably been observed in [15]. However, in [16] two additional 
states, lying at rather low energies, have been observed which do not seem to correspond to the 
calculated levels. It has to be mentioned that neither a shell model calculation, performed within 
an harmonic oscillator basis, provides satisfactory results in this case. Yet, we found that this 
shell model calculation works better than one would assume given the unstable character of the 
states involved. In four-particle system we found that there is no bound or antibound state in 
13 Li, which agrees with the experiment that the 13 Li is unbound. 
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